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Abstract The Erdos-Renyi classical random graph is characterized by a fixed linking 
probability for all pairs of vertices. Here, this concept is generalized by drawing the linking 
probability from a certain distribution. Such a procedure is found to lead to a static complex 
network with an arbitrary connectivity distribution. In particular, a scale-free network with the 
hierarchical organization is constructed without assuming any knowledge about the global 
linking structure, in contrast to the preferential attachment rule for a growing network. The 
hierarchical and mixing properties of the static scale-free network thus constructed are studied. 
The present approach establishes a bridge between a scalar characterization of individual 
vertices and topology of an emerging complex network. The result may offer a clue for 
understanding the origin of a few abundance of connectivity distributions in a wide variety of 
static real-world networks. 
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In recent years, a lot of efforts have been made to characterize the structures of real-world 
networks observed in physical, biological, economical, and man-made systems [1-3]. An 
intriguing feature is that most of these networks turn out to belong to one of the following 
classes: small-world [4], scale-free [5], exponential, or random [6]. Note that although all of 
these types may involve intrinsic randomness, referred to as the "random graph" is the 
classical one introduced by Erdos and Renyi [6]. 

While the preferential attachment rule [2,5] can nicely explain how a growing network 
self-organizes to scale-free topology, there exist a number of real scale-free networks, which 
do not follow such a mechanism. Some examples are given in Ref. [7]. Of particular, interest 
is the static scale-free network, in which the total number of vertices does not change in time, 
i.e., no growth. It is therefore necessary in a coherent setting to develop a theory of static 
complex networks and to examine their connectivity distributions, clustering properties, and 
levels of hierarchy [4,8]. In such a direction, some fruitful discussions have been developed in 
the literature. In Refs. [9-11], a static model of varying vertex fitness has been introduced and 
analyzed. In analogy with the superstatistics formulation [12,13] of nonextensive statistical 
mechanics [14,15], a general analytic result has been presented in Ref. [16] for the relation 
between the linking probability distribution for a fluctuating random graph (i.e., varying 
fitness) and a complex network with the connectivity distribution of an arbitrary form. 
Another approach to static complex networks has been proposed in Refs. [17,18], in which a 
constant-size scale-free network is made by imposing the balance between creation and 
merger of vertices. 

In this paper, we study the seemingly simplest modification of the Erdos-Renyi classical 
theory of random graphs. The classical random graph is constructed by using the linking 
probability which is kept fixed for all pairs of vertices. Here, we start from a statistical 
characterization of the linking probabilities between vertices, and arrive at static complex 
network with various types of connectivity, including scale-free topology. We shall not 
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employ the ensemble picture in Ref. [16] but introduce an algorithm, which enables one to 
produce a single realization of each network. We are motivated by the observation of the vast 
abundance of similar connectivity distributions in a wide variety of static real-world networks, 
in which vertices may be individuals, animals, chemicals, companies and so on. Often, 
vertices are characterized by a linking probability defining the way of connecting them by 
edges. A linking probability can be nonuniform for all pairs of vertices, in general, unlike in 
the classical theory of random graphs. It conditions characteristics and states of vertices, and. 
introduces a nontrivial correlation property. We may analyze the hierarchical and mixing 
properties of the static complex networks thus constructed. 

The present approach may offer a possibility to find a common "driving force" leading to 
the fact that only a few classes of networks are observed in the real world. This driving force 
would be featured by a linking probability distribution common in various systems. For 
example, many social and economic networks explicitly depend on wealth of individuals, 
firms, banks, and so on. Suppose that the wealth distribution, which obeys a power law 
common in industriarized countries [19,20], is associated with such a linking probability 
distribution, and assume that, other than this probability, linking is entirely random. Then, the 
resulting network will have a connectivity distribution, which is uniquely determined by the 
linking probability distribution [16]. A common driving force could thus provide an 
understanding of why real static networks tend to belong to only a few "universality classes". 

Let us start our discussion with the classical theory of random graphs presented in the 
milestone papers of Erdos and Renyi [6]. There are at least two known ways to construct 
random graphs with N vertices. The first one is to connect any two of all possible 
A'^( A'^ - 1) / 2 pairs of vertices with probability p s (0, 1) and not to connect with probability 
I- p. Then, connectivity is described by the binomial distribution. 



p" (I- p)'^ ^ which in the large A'^ limit becomes the Poissonian 



distribution 



-X 
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where A is the average value of connectivity, X=<k>= p(N-l), and is to be kept fixed. 
This distribution has finite moments and a characteristic scale, A. The second one is to chose 
L edges at random and connect their end vertices to any two vertices. In this case, A 
corresponding to that in the first construction is given by <k>=2L I N . The latter 
construction allows the presence of tadpoles (i.e., self-interaction) and melons, in general. In 
the large limit with fixed A, however, these two methods are statistically equivalent, and 
the resulting connectivity distributions are Poissonian as in Eq. (1). 

Besides the connectivity distribution, there are several other quantities characterizing 
networks. One such is the clustering coefficient, c{k), defined by 



where A:, is connectivity of the /th vertex and P{k) the connectivity distribution of the 

network (as a simple graph) under consideration. The quantity is related to the total 

number of the edges connecting all the nearest neighbors of the ith vertex. Using the 
adjacency matrix A = (a,^) [that is, a,^ = 1 (a,^ = 0) if the rth and j\h vertices are directly 

connected (unconnected), and a^, = because of the absence of loops], it can be expressed as 




(2) 



c, = 



k,{k,-\y 



(3) 



4 



(4) 



For the Erdos-Renyi classical random graphs as well as the scale-free networks generated by 
growth with the preferential attachment rule (i.e., the Barabasi- Albert networks [5]), c(k) 
does not depend on k, whereas the so-called hierarchical networks [8] have c(k) decaying as 
a power law with respect to k. The global clustering coefficient C is defined by the average of 
c- over all vertices [4]. The small-world networks have the values of C, which are much 

larger than those of the Erdos-Renyi classical random graphs [4]. 

Another useful quantity is the correlation coefficient between connectivities of nearest- 
neighbor vertices [3,21]: 

k„„(k) = Y,k' P(k'\k), (5) 

k' 

If k„„(k) increases (decreases) with respect to k, the mixing property [3,21,22] is called 

assortative (disassortative): that is, vertices with large values of connectivity tend to have a 
majority of neighboring vertices with large values (small values) of connectivity. It is known 
[22] that the networks of scientific coauthorships and actor collaborations have assortative 
mixing, whereas the Internet, the World Wide Web, protein architectures in cells, and neural 
networks possess disassortative mixing. Very importantly, neither the Erdos-Renyi classical 
random graphs nor the Barabasi-Albert scale-free networks show these properties, that is, 
their k^^(kys are actually independent of k. 

Once again, we notice that it is a crucial point in the Erdos-Renyi classical random graphs 
that each vertex has the same linking probability to get connected with one of the A^-1 
potential neighbors. Now, suppose that the linking probability is not fixed but can vary for 
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each pair of vertices. In this case, p, or equivalently A in Eq. (1), may be a random variable 
obeying a certain probability distribution IT (A), which is sometimes referred to 
(inappropriately) as the "hidden variable distribution". Then, P^^{k) in Eq. (1) has to be 

regarded as a conditional probability distribution. Then, the correct connectivity distribution 
P{k) is given by the marginal distribution of the joint distribution, n(A) P^Ak) [15,23]: 



This has a nice parallelism with the superstatistics approach to nonextensive statistical 
mechanics [12-15]. Conventional statistical mechanics characterized by the Boltzmann factor 
is a theory for simple systems in chaos, whereas nonextensive statistical mechanics is 
designed for complex systems at the edge of chaos. Superstatistics introduces significant 
fluctuations of the local inverse temperature jS, so that the temperature in the Boltzmann 
factor is stochastic and obeys a certain distribution function (commonly the chi square 
distribution). The resulting superposed distribution turns out to be that of nonextensive 
statistical mechanics. Thus, we see the following correspondence relations: (Boltzmann-Gibbs 
statistical mechanics) (random graphs), (nonextensive statistical mechanics) (complex 
networks), and (fluctuating fi) (fluctuating A). (For another consideration of the relation 

between complex networks and nonextensive statistical mechanics, see Ref. [24].) 

Now, the problem is to find IT (A) for a given P{k). An analytic method for solving this 
has recently been presented in Ref. [16]. In particular, consider a static scale-free network 
characterized by the power-law connectivity distribution 




(6) 
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with A ' = ^(7,^0), ^oG(0, 1), 7>1, and Hurwitz's generalized zeta function ^(s,a). 
Then, the exact fonn of IT (A) is given by [16] 



where r(z) is Euler's gamma function. 

Now, let us describe our algorithm for the above general construction not in employing the 
ensemble picture but in aiming to construct a complex network in each realization, in contrast 
to the standpoint in Ref. [16]. To be manifest, recall the following procedure of generating the 
Erdos-Renyi classical random graph: (i) fix the total number of vertices, A^, (ii) fix the linking 
probability, p, for all vertices, that is, fix X=<k>= p{N -V) d& well as the total number of 
edges, L = pN{N -\) 1 2 = XN 1 2, and then, (iii) pick up the /th and 7th vertices with the 
same probability p, connect them, and repeat this L times. We now generalize this algorithm 
as follows: (I) fix the total number of vertices, N, (II) fix the linking probability distribution 
for all vertices, that is fix A,- = p^{N -I) for the iih vertex as well as the total number of 

edges, L=< p.> N(N -I) / 2=< Xf> N / 2 with the brackets being the averages with 
respect to the linking probability distribution [the distributions, Il(p) and FT (A), being 
different from each other only by a constant scale factor of the variables, A'^-1], and then 
(HI) pick up the ith and 7th vertices with the probabilities p- and Pj, respectively, connect 

them, and repeat this L times. 

To realize the algorithm (I)-(III), we proceed in the following way. We assign the random 
variable p- drawn from the distribution Il(p) to all vertices. Then, we impose the 

normalization condition, 2^._^ Pi = 1, eliminating any spurious influence of the constant 




(8) 
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scale factor, A'^-l. When connecting vertices, the probability that the ith vertex can be 
connected to one of two ends of the edge is . Here, as an example, we allow a vertex to be 

connected to itself twice, leading to tadpoles. This is analogous to the second method of 
constructing the Erdos-Renyi classical random graph with tadpoles and melons, mentioned 
just after Eq. (1). We have numerically ascertained that elimination of tadpoles does not affect 
the results reported below for the network sizes A'^ = 1000 ~ 10000. 

In Fig. 1, we present three different forms of 11 (A): 11 (A) = 5(X-ri) with constant rj, 

n(A,) = ae'"^ (a> 0), and n(A) = 11^^ (A) in Eq. (8), which lead to the Poissonian (Erdos- 
Renyi classical random graph), exponential (exponential network), and power-law (scale-free 
network) connectivity distributions, respectively. As can be seen from Eq. (8), 11 (A) decays 

as a power law with the same exponent as that of the connectivity distribution in Eq. (7). 

Among the above three types, most interesting may be the scale-free network. In Fig. 2, we 
show the plots of the clustering coefficient of the scale-free network for the three different 
values of the exponent: 7 = 1.5, 2.0, 2.5. c(k) has been calculated for each realization of the 
network, using the algorithm explained earlier and analyzing the corresponding adjacency 
matrix. 

There is a technical point in calculating c(k) when a given realization does not contain 
vertices with some values of connectivity. To treat such a case, we have examined the 
following two averaging methods. In the first one, averaging is carried out only for the values 
of connectivity present in the realization. The result obtained by this method is shown in Fig. 
2 (a). The second method is that the clustering coefficient is set equal to zero for values of 
connectivity absent in the realization and then averaging is performed. The result based on 
this method is presented in Fig. 2 (b). In both cases, decay of c(k) with respect to k is 
observed, manifesting hierarchical organization of the present static scale-free networks. It is 
quite remarkable that 11^^ (A) in Eq. (8) naturally leads to such a structure, without employing 
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any artificial treatments such as thresholding (see, e.g., Ref. [10]). 

We have also ascertained that the Erdos-Renyi classical random graph corresponding to the 
top one in Fig. 1 has c(k) that does not depend on k, reproducing the well known fact that 
that a classical random graph does not possess the hierarchical structure. However, we are not 
able to arrive at definite conclusions about the hierarchical structure of the exponential 
networks corresponding to the middle one in Fig. 1, since it turned out to be numerically very 
hard for our current computational power to obtain reasonable statistics for large values of 
connectivity, unfortunately. 

Finally, in Fig. 3, we present the result about the mixing property of the static scale-free 
network. Corresponding to the power-law connectivity distributions in Fig. 1 (c), the 
correlation coefficient between connectivities of nearest-neighbor vertices, (k) , in Eq. (5) 

has been calculated for different values of the exponent of the connectivity distribution: 
7 = 1.5, 2.0, 2.5. The decay of correlation is clearly appreciated. Consequently, we find that 
the static scale-free network constructed here possesses disassortative mixing. 

In conclusion, we have studied generalizations of the Erdos-Renyi classical theory of 
random graphs to allow variable linking probabilities for pairs of vertices. In particular, we 
have discussed in detail the case of static scale-free networks and have seen how hierarchical 
organization and disassortative mixing are realized. We think that special cases of the 
algorithm presented here might have partially been considered in Refs. [10,11]. 

The present study is based on a recent work in Ref. [16], in which it is analytically shown 

how to construct, from the classical random graph, a static complex network with an arbitrary 

connectivity distribution. It is of interest to see that there exists a bridge between the global 

topological properties (e.g., connectivity and clustering) and the local individual ones (i.e., 

linking probabilities). This may potentially help to explain why diverse real- world networks 

have common features. Such a view is particularly attractive if the fact that the approach does 

not rely on global information about the entire network (such as the preferential attachment 
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rule or its variants) is taken into account. 
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Figure captions 



Fig. 1 The connectivity distributions for three different linking probabilities with 
A'^ = 10000 vertices. The results are obtained by averaging over 10 identical network 
realizations. Top: the Poissonian connectivity distributions of the Erdos-Renyi classical 
random graphs corresponding to the Dirac delta distribution for IT (A). Middle: the 
exponential connectivity distributions corresponding to the exponential distributions for 
n(A). Bottom: the power-law connectivity distributions of the scale-free networks 
corresponding to (A) in Eq. (8). The solid lines represent the analytical results obtained in 
Ref. [16]. 
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Fig. 2 The log-log plots of the dependencies of the clustering coefficient of the static scale- 
free networks on connectivity for three typical values of the exponent of the power-law 
connectivity distribution obtained by two different methods: (a) averaging is carried out only 
over connectivity present in each network realization, and (b) the clustering coefficient of 
connectivity absent in each network realization is set equal to zero and then averaging is 
carried out. From top to bottom in both (a) and (b): 7 = 1.5, 2.0, 2.5. Averaging is performed 
over 100 independent realizations. Both (a) and (b) yield non-trivial hierarchical organization. 
Here, averaging does not mean the ensemble picture: each configuration has these behaviors 
(i.e., averaging is nothing but for obtaining statistically rather smooth curves). 

Fig. 3 The log-log plots of the correlation coefficient between connectivities of nearest- 
neighbor vertices, k^(k), with respect to connectivity for three typical values of the 

exponent of the power-law connectivity distribution. From the top curve to the bottom 
one: 7 = 1.5, 2.0, 2.5. Averaging is performed over 100 independent realizations (which are the 
same configurations as in Fig. 2) using the linking probability 11^,^(1) in Eq. (8), for 

networks with A'^ = 1000 vertices. 
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